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ABSTRACT 



We combine the CfA Redshift Survey (CfA2) and the Southern Sky Redshift Sur- 
vey (SSRS2) to estimate the pairwise velocity dispersion of galaxies (J12 on a scale of 
~ l/i~^Mpc. Both surveys are complete to an apparent magnitude limit B{0) — 15.5. 
Our sample includes 12,812 galaxies distributed in a volume 1.8 x 10'^ /i~^Mpc^. We con- 
clude: 

1) The pairwise velocity dispersion of galaxies in the combined CfA2+SSRS2 redshift 
survey is (J12 = 540 km s"-*- it 180 km s"-*-. Both the estimate and the variance of (J12 
significantly exceed the canonical values (J12 = 340 it 40 measured by Davis & Peebles 
(1983) using CfAl. 

2) We derive the uncertainty in (J12 from the variation among subsamples with volumes 
on the order of 7 x 10^ Mpc^. This variation is nearly an order of magnitude larger 
than the formal error, 36 kms~^, derived using least-squares fits to the CfA2+SSRS2 
correlation function. This variation among samples is consistent with the conclusions of 
Mo et al (1993) for a number of smaller surveys and with the analysis of CfAl by Zurek 
et al. (1994). 

3) When we remove Abell clusters with R > 1 from our sample, the pairwise velocity 
dispersion of the remaining galaxies drops to 295 km s"-*- it 99kms~-'-. 

4) The dominant source of variance in a 12 is the shot noise contributed by dense 
virialized systems. Because (J12 is pair- weighted, the statistic is sensitive to the few richest 
systems in the volume. This sensitivity has two consequences. First, a 12 is biased low 
in small volumes, where the number of clusters is small. Second, we can estimate the 
variance in a 12 as a function of survey volume from the distribution of cluster and group 
velocity dispersions n{a). For either a COBE-normalized CDM universe or for the observed 
distribution of Abell cluster velocity dispersions, the volume required for (J12 to converge 
(^0-12/0-12 < 0.1) is ~ 5 X lO^/i-3 Mpc^, larger than the volume of CfA2+SSRS2. 

5) The distribution of pairwise velocities is consistent with an isotropic exponential 
with velocity dispersion independent of scale. The inferred szra^f/e-galaxy velocity distribu- 
tion function is incompatible with an isotropic exponential. Thus the observed kinematics 
of galaxies differ from those measured in the hydrodynamic simulations of Cen & Ostriker 
(1993). On the other hand, our observations appear to be consistent with the velocity 
distribution function measured by Zurek et al. (1994) using coUisionless simulations with 
much higher resolution than Cen & Ostriker (1993). The large dynamic range in the sim- 
ulations by Zurek et al. (1994) affords a more accurate treatment of galaxy interactions, 
which dramatically alter the dynamical evolution the galaxy distribution (Couchman & 
Carlberg 1992, Zurek et al. 1994). The agreement between the observed velocity distribu- 
tion and the one predicted by these high-resolution simulations may be another clue that 
mergers play an important role in the evolution of galaxies and large-scale structure. 
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1 INTRODUCTION 



Galaxy motions offer a glimpse of dark matter in the universe. Although individual 
peculiar velocities are notoriously difficult to measure directly, relative velocities of galaxy 
pairs can be measured from redshift surveys alone. Redshifts combine the cosmological 
component Hor with the line-of-sight peculiar velocity. Because the distribution of galax- 
ies on the plane of the sky is unaffected by peculiar motions, correlation functions (or, 
equivalently, power spectra) which are isotropic in real space are anisotropic in redshift 
space {e.g., Sargent & Turner 1977, Peebles 1979, Peebles 1980, Kaiser 1987, Hamilton 
1992). The degree of anisotropy in redshift space measures the low-order moments of the 
peculiar velocity distribution. 

These moments play an important role in the evolution of large-scale structure. The 
first moment of the pairwise velocity distribution, Ui2, governs the growth of the spatial 
correlation function through the conservation of particle pairs (Davis & Peebles 1977, 
Peebles 1980, Fisher et al. 1994b). If galaxies trace mass and density fiuctuations are in 
the linear regime, then Vi2 is proportional to il^-^. Fisher et al. (1994b) demonstrate a 
technique to derive Vi2 from redshift surveys using N-body simulations, and then measure 
Vi2 from the IRAS 1.2 Jy survey. 

The second moment of the pairwise velocity distribution, (J12, measures the kinetic 
energy of random motions in the galaxy distribution; in equilibrium, this dispersion ex- 
actly balances the gravitational potential. (Geller & Peebles 1973, Peebles 1976, Davis & 
Peebles 1977). Peebles (1976) measured the anisotropy of the redshift-space correlation 
function and derived (J12 from the Reference Catalog of Bright Galaxies (deVaucouleurs & 
deVaucouleurs 1964). Davis & Peebles (1983, hereafter DP83) extended this analysis to a 
much larger redshift survey, CfAl (Huchra et al. 1983). DP83 found an exponential pair- 
wise velocity distribution with (J12 = 340 it 40kms~^. This observation has since become 
one of the strongest constraints on models for the evolution of large-scale structure, and 
has remained the standard against which N-body simulations are judged. 

The favored paradigm for the origin of structure follows the gravitational collapse of 
dense regions from initially small perturbations to the large, non-linear structures we ob- 
serve today. Among the many possible constituents of the initial density field, cold dark 
matter (CDM) stands as the theoretical standard of comparison. In its standard form, 
CDM presumes a fiat universe; the only free parameter in the model is the normalization 
of the density power spectrum, usually cast as the rms density fiuctuation in S Mpc 
spheres, erg. Both the quadrupole anisotropy in the CMB (Wright et al. 1992) and the 
distribution of optically selected galaxies (Vogeley et al. 1994, Loveday et al. 1993) suggest 
that (jg ~ 1. Using this normalization, N-body simulations of a fiat CDM universe yield 
(712 ~ 1000 km s~^, much larger than the values obtained by DP83. Open {ilh — 0.2) mod- 
els predict (J12 ~ 500kms~^, thus reducing but not eliminating the discrepancy (Davis et 
al. 1985). Under the simple assumption that galaxy velocities trace the mass distribution, 
the apparent quiescence of the observed velocity field strongly constrains CDM models. 

A fair comparison between theory and observation requires an accurate prescription 
for galaxy formation. The possibility remains that the galaxies we observe paint a biased 
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picture of the matter distribution as a whole. If, for example, galaxies form only in dense 
regions, then galaxies naturally cluster more than the underlying density field (Kaiser 1984, 
Bardeen et al 1986). This density bias lowers the normalization of the mass power spec- 
trum corresponding to the observed galaxy power spectrum, and thus lowers the expected 
peculiar velocities. In simulations of a biased CDM universe, Davis et al (1985) calculate 
(712(0. Mpc) ~ 700 kms-^ 

The COBE normalization is more restrictive. Anisotropies in the microwave back- 
ground trace fiuctuations in the potential itself; thus the COBE measurement removes 
the freedom in erg allowed by density bias alone. However, Carlberg (1991) suggests that 
galaxy velocities may be biased as well. Carlberg (1991,1994), Carlberg & Couchman 
(1989), Evrard et al. (1992), Katz et al. (1992), and Cen k Ostriker (1993) describe a 
number of processes which rob energy from galaxy tracers either through purely gravi- 
tational interactions or through dissipation in the baryonic component. On the whole, 
these biases reduce the simulated (J12 by up to a factor of ~ 2 on scales smaller than a 
megaparsec, not quite enough to bring the models into agreement with DP83. 

The range of theoretical possibilities prompts a re-examination of the observed (J12. 
Table 1 summarizes earlier work. The variations in the observed a 12 are large. Although 
the largest surveys to date, CfAl (Huchra et al 1983) and IRAS12 (Fisher et al 1994b), are 
consistent with each other to within the quoted errors. Mo et al (1993) find that estimates 
of a 12 from surveys of this size are extremely sensitive to the treatment of the few rich 
clusters they contain. Zurek et al (1994) reach similar conclusions in their comparison 
between CfAl and a series of large, high-resolution simulations. Zurek et al (1994) find 
substantial fiuctuations among simulated volumes on the scale of CfAl. Furthermore, 
both Mo et al (1993) and Zurek et al (1994) claim that CfAl is, in fact, consistent with 
simulations if the redshifts are not corrected for infall into the Virgo cluster. 

Here, we use a sample of 12,812 galaxies drawn from the CfA Redshift Survey (Huchra 
et al 1995, CfA2 hereafter) and the Second Southern Sky Redshift Survey (da Costa et 
al 1994, SSRS2 hereafter) to measure (J12 and to examine the variance of (J12 in volumes 
of ~ 10'^/i~^Mpc^. In §2, we describe the data. In §3, we review the method we use to 
extract (J12 from the anisotropy of the correlation function. We test the method in §4 using 
N-body and Monte-Carlo simulations of the galaxy distribution. Section 5 summarizes the 
measurements of (J12 using the distribution functions discussed in §3. In §6, we discuss 
the implications for cosmological models and examine the convergence of (J12 with survey 
volume by relating (J12 to the distribution of group and cluster velocity dispersions. We 
conclude in §7. 
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2 DATA 



The extension of the CfA Redshift Survey (deLapparent et al 1986, Geller & Huchra 
1989, Huchra et al 1994, hereafter CfA2) to B{0) < 15.5 now includes approximately 
11,000 galaxies. We analyze two subsamples, one in the north Galactic cap (CfA2N) and 
one in the south Galactic cap (CfA2S). CfA2N contains 6,480 galaxies over the region 
S'^ < a < 17^ 8.5' < 6 < 44. 5^ CfA2S includes 4095 galaxies over the region 20^^ < 
a < 4^ and —2.5'' < 6 < 42°. Because the survey was incomplete when we began this 
analysis, we do not include the northernmost slice of CfA2S included in Huchra et al 1994. 
Because CfA2S suffers heavy extinction along its boundaries, we exclude the following 
regions from our analysis (cf. Vogeley et al 1994): 21'^ < a < 4^;21^ < a < 2^ and 
jyii > -25^2^^ < a < 3^^ and b^^ > -45% where b^^ is the galactic latitude. The 
remaining CfA2S sample covers 0.77 steradians and includes 2741 galaxies. 

We draw CfA2 galaxies from the Zwicky catalog. Although better photometry exists 
for a small fraction of these galaxies, we use the original Zwicky magnitudes for all galaxies. 
A small fraction of the galaxies in our sample are multiple systems for which Zwicky et 
al (1961-1968) estimated only a combined magnitude. For these cases, we estimate the 
relative contributions of each component by eye, and eliminate galaxies which fall below 
the magnitude limit. 

We correct heliocentric redshifts for the solar motion with respect to the centroid of 
the Local Group: Av — 300 sin (/) cos (6) km s"-*-. We do not correct for infall into the Virgo 
cluster. We exclude galaxies with redshifts smaller than 500 kms~^ and larger than 12,000 
kms"-*-. Beyond 12,000 kms"-*-, fewer than 1% of the galaxies visible at 500 kms"-*- are 
brighter than the magnitude limit given the luminosity function for CfA2N (Marzke et al 
1994). The mean uncertainty in the measured redshifts is 35 kms~^. 

The Second Southern Sky Redshift Survey (da Costa et al 1994, hereafter SSRS2) 
includes 3591 galaxies over the region —40'' < 6 < —2.5°, b^^ < —40°. Alonso et al. 
(1993a, 1993b) selected the sample from a combination of two machine-generated two- 
dimensional catalogs: the Automatic Plate Measuring (APM) Galaxy Survey (Maddox et 
al. 1990) and the non-stellar object list of the STScI Guide Star Catalog (GSC)(Lasker et 
al. 1990, and references therein). Alonso et al. (1993a) calibrated the magnitude scale using 
overlapping galaxies in the Surface-Photometry Catalogue of the ESO-Uppsala Galaxies 
(Lauberts & Valentijn 1989). The ESO catalog also provided the bright galaxies missed by 
the automatic identification algorithms; these galaxies comprise 10% of the sample. For 
the SSRS2, the mean uncertainty in the redshifts is ~ 50kms~^. 

These samples demonstrate rich structure in the nearby universe (Geller & Huchra 
1989, da Costa et al. 1994). CfA2N contains the Great Wall, which includes the Coma 
cluster at ~ 7, 000 km s~^. This sample also includes the Virgo cluster nearby and several 
large voids in between. The most striking feature of CfA2S is the dense wall connecting 
the Perseus and Pisces clusters. A similar wall appears in the SSRS2 running diagonally 
across the survey. 

Table 2 summarizes the basic features of these samples. The combined CfA2+SSRS2 
sample includes 12,812 galaxies in an effective volume of 1.8 x 10'' Mpc^. 
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3 METHOD 



In this section, we outline the method used to extract the pairwise velocity dispersion, 
(712, from redshift survey data. The technique was pioneered by Geller & Peebles (1973) and 
honed by Peebles (1976,1979), Davis et al (1978) and Davis k Peebles (1983). Briefly, the 
method requires three steps: 1) establishing the variation of the redshift-space correlation 
function (CF, hereafter) with angle to the line of sight, 2) measuring the spatial CF from 
the projected separations and mean redshifts of pairs and 3) modeling the redshift-space 
CF as a convolution of the spatial CF with the distribution of relative peculiar velocities 
for pairs of galaxies. 



3.1 Definition and Computation of^{rp,Ti) 

To isolate distortions in the CF caused by peculiar velocities, we divide the pair sepa- 
ration in redshift space, s, into two components, one along and the other perpendicular to 
the line of sight to the pair (Davis, Geller & Huchra 1978). The line of sight, 1, bisects the 
pair: 1 = ^(si + S2) where Sj- = czi^i/ Hq, c is the speed of light and the redshift z <^ 1. 
Assuming that the angular separation of the pair is small, the components of s = si — S2 
are 

si 

TT = — — - 



p 



s|2-7r2 (1) 



Throughout this paper, we express the line-of-sight component vr in units of distance. For 
pairs with small angular separations, vr measures a combination of the Hubble separation, 
ri2, and the component of the pairwise velocity ,vi2, along ri2: tt |ri2| + (vi2 • ri2/|ri2|). 
Following DP83, we restrict our analysis to pairs separated by angles smaller than 9max — 
50''. This choice of 9max balances the need for densely sampled data against the fact that 
widely separated pairs overestimate V12 in the plane-parallel approximation (1). 

The two-dimensional CF ,C(rp,7r), is deflned by the probability that any two galaxies 
in the survey are separated by Vp and vr (Peebles 1980): 

6P ^ n[l + ^{rp,Tr)]6Tr6A (2) 

where 6A — 'lurpdrp. To measure C(rp, vr), we compute the distribution of pair separations 
in the data and in a Poisson realization of the data with the same radial and angular 
selection criteria (Peebles 1980): 

DR{rp,Ti)nD 
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where DD and DR are the number of data/data and data/random pairs, respectively, with 
separations Vp and vr. We compute the mean densities of the real and random samples, n/j 
and ur, using the minimum- variance estimator derived by Davis & Huchra (1982). 

In order to use the estimator in Equation 3, we need to choose appropriate bins in Vp 
and TT. Linear binning in both dimensions provides a clear picture of various contributions 
to the redshift-space anisotropy (cf. Fisher et al 1994b). On the other hand, our prior 
knowledge of the form of the CF in Vp and vr suggest a mixture of logarithmic bins in Vp and 
linear bins in vr would be more appropriate (DP83). For our analysis, we use logarithmic 
bins in Vp and linear bins in it. In order to facilitate comparisons with other surveys, we 
display the CFs using linear bins in both dimensions; again, we use this binning for display 
purposes only. We compute the correlation function on a 30x30 grid with 1 Mpc bins in 
TT and O.ldex bins in Vp. The 30 bins in Vp cover the range 0.1 Mpc to 100 Mpc. To fit 
models to C(rp,7r), we use coarser binning in Vp to increase the signal-to-noise: each bin 
represents a factor of 2 in Vp starting at Tp = 0.1 Mpc. 

To construct the random samples, we use the luminosity functions derived by Marzke 
et al (1994) for the CfA survey and by da Costa et al (1994) for the SSRS2. Because 
we lack morphologies for most galaxies in the sample, we use the general luminosity func- 
tion averaged over all types. The luminosity functions for CfA2N, CfA2S and SSRS2 
differ significantly. In order to remain self-consistent, we construct the random sample 
for CfA2+SSRS2 using the individual luminosity functions rather than using a mean av- 
eraged over all samples. Although these variations in the luminosity function complicate 
the interpretation of large-scale features in C(rp,7r), they do not affect our analysis of the 
small-scale velocity field. 

To compute the pair sums DD and DR, galaxies may be counted in any number of 
ways as long as the weight assigned to each galaxy does not depend explicitly on the local 
density (Hamilton 1993). In a "fair" sample, almost by definition, any weighting yields 
an unbiased estimate. The simplest approach is to count each galaxy equally. Because 
the density of galaxies in a magnitude-limited sample drops rapidly with distance, this 
weighting favors nearby galaxies. Weighting by the integral of the luminosity function 
to account for the magnitude limit gives equal weight to each volume of space (Davis & 
Huchra 1982). In the minimum- variance (MV, hereafter) estimate of ^(rp,7r) for a fair 
sample of galaxies (Efstathiou et al 1988), each pair carries a weight wiW2 where 



1 + 4:Trn(f){si)J3{sij) 

Here, Si is the redshift of galaxy i, Sij is the pair separation in redshift space, (f) is the 
selection function and Js is the volume integral of the spatial CF (Peebles 1980). MV 
weighting balances the advantage of equal- volume weighting against the shot-noise caused 
by the few close pairs at large redshift (Efstathiou et al 1988, Loveday et al 1992, Park 
et al. 1994,Fisher et al. 1992). 

The scale of the CfA2+SSRS2 redshift survey is no larger than the scale of the largest 
apparent structures. In these samples, different weighting schemes yield systematically 
different estimates of the CF, depending on the variation of clustering with distance and 
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with galaxy luminosity (deLapparent et al 1988, Park et al 1994). We show in §5 that 
variations in C(rp,7r) derived with different weighting schemes are not negligible. Because 
MV weighting ensures the most efficient use of the data, we present our final numerical 
results using this weighting exclusively. 

Correlation statistics of samples containing structure on the scale of the survey are 
particularly sensitive to uncertainty in the mean density of the universe (deLapparent et 
al 1988, Kaiser 1987, Hamilton 1993). Because the mean is most often defined internally 
as the number of sample galaxies divided by the sample volume, the volume integral of 
the CF vanishes artificially on the scale of the survey. On large scales, where the CF is 
small, errors in the CF are proportional to the error in the mean density, 6 (deLapparent 
et al 1988). Hamilton (1993) suggests an alternative to the standard estimator (Equation 
3) which is less affected by uncertainty in the mean density: 



i+Mrp,-) = ^^^;= (5) 



{ DD){R R) 

Because the error in is proportional to d'^ , this estimator performs better than Equation 
3 when 6 <^1. In the samples we analyze here, 6 is not much smaller than unity (Marzke 
et al. 1994, da Costa et al. 1994), and ~ For our purposes, the differences between 
the estimators in Equations 3 and 5 are negligible. Fisher et al. (1994b) reach similar 
conclusions in their analysis of the IRAS 1.2Jy survey. 



3.2 Model for ^{rp,Tr) 

^{rp, it) may be expressed as a convolution of the spatial CF, C(r), with the distribution 
of pair wise peculiar velocities (Peebles 1980): 

l + e(rp,7r)= f[l + ar)]P{^i2\ri2)d^^i2 (6) 

If the pairwise peculiar velocity distribution varies slowly with pair separation and if there 
is no preferred direction in the velocity field, then C(rp,7r) may be expressed in terms of 
the distribution of line-of-sight velocities alone: 

l + e(rp,7r) = J[l + ar)]P{vjos\ri2)dvjos (7) 

If we similarly decompose the spatial separation ri2 into components (rp, y) perpendicular 
to and along the line of sight, then r^2 — ''r v"^ , ^los — — y and 

l + e(rp,7r) = J[l + ^(^r^^ + y2)]p((^-y)\r,,)dy (8) 

(Peebles 1980). 

This expression is particularly simple when both members of the pair share the same 
line of sight, i.e. Tp = 0. In this case, 

1 + arp = 0, tt) = y'[i + avMi^ - y)\y)dy (9) 
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Thus if we know the spatial CF and we measure ^{rp, it), we can probe the peculiar velocity 
distribution. Specifically, we measure the distribution function of a single component of 
the pairwise velocity, P{vi2)- We will see later that the cut in C(rp,7r) at Tp = best 
discriminates among models for P{vi2)- 



3.3 Distortions m ^{rp,Ti) from Peculiar Velocities 

Two limiting cases illuminate the distortions in C(rp,7r) caused by peculiar velocities. 
On scales where the width, (J12, of the pairwise velocity distribution exceeds the Hubble 
separation (for example in relaxed clusters and groups), pairs that are close in space are 
on average more widely separated in redshift space. Thus the slope of the redshift-space 
CF is shallower along the line of sight than perpendicular to it; the contours of C(rp,7r) 
are extended along the line of sight. This distortion is the familiar "redshift finger", a 
distinctive feature of galaxy redshift surveys. The length of the fingers is a standard 
measure of the second moment of the velocity distribution function. 

Another distortion is caused by the mean motion of galaxies toward one another. 
On scales where the mean infall velocity is smaller than the Hubble separation, infall 
compresses C(rp,7r) along the line of sight (Peebles 1980). Because density fiuctuations 
on these scales are small, linear perturbation analysis yields an expression relating the 
amplitude of this distortion to the mean mass density of the universe Q, (Peebles 1980, 
Kaiser 1987, Hamilton 1993, Fisher et al. 1994b). Given a model for the bias between 
galaxies and mass, the mean infall velocity provides a direct route to il. Unfortunately, 
coherent structure on the scale of the survey masquerades as large-scale streaming (Fisher 
et al. 1994b), and the true infall velocity is nearly impossible to disentangle. 



3.4 Spatial Correlation Function 

Our ability to explore the pairwise velocity distribution through Equation 8 rests on 
our knowledge of the true distribution of galaxies in space. If we measure C(rp,7r) on 
large enough scales, we can deproject ^(r) from C(rp,7r) without modeling the velocity 
distribution function at all (DP83). The projection of ^(rp,7r) onto the Tp axis, w{rp), 
depends only on the spatial CF: 

J 

= / dy^{^r^ + y^) r-p < ir^ax (10) 

In practice, we measure C(rp,7r) to a scale iTmax where the signal-to-noise approaches 
unity. The projection in Equation 10 is independent of the velocity distribution only if 
T^max ^ Ci2- Given w{rp), we invert Equation 10 for ^(r) (DP83): 

a.) = -iri5^ (11) 

TV Jr , / f ^ — f ^ 
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For a power law spatial CF ^(r) = {r/roy , the projected CF is w{rp) — Arp~^^ with 



^,r(l/2)F[(T-l)/2] 



r(7/2) ^ ' 



3.5 Velocity Distribution Function 



In linear theory, the peculiar velocity field of individual galaxies is simply related to 
the density field (Peebles 1980): 

Vl r = ^ / ^ ^ d^r' 13 

If the one-point density distribution is Gaussian, then the velocity distribution function 
is Gaussian and is completely specified by the power spectrum of density fiuctuations 
(Vittorio ct al 1989, Scherrer 1992): 

l-CO 

al = Hln^ '^ / 4TrP(k)dk (14) 
Jo 



On the small scales we observe here, non-linear evolution of the phase-space distribution 
erases much of the history written in the initial galaxy orbits. In this regime, the velocity 
distribution is difficult to predict. Our knowledge of the small-scale velocity distribution 
comes primarily from N-body simulations and is further guided by observations. Peebles 
(1976) calculated C(rp, vr) using the Reference Catalogue of Bright Galaxies (deVaucouleurs 
and deVaucouleurs 1964) and observed that an exponential distribution of pairwise line- 
of-sight velocities best reproduced the data: 

F«2) = exp(-y2<2/cTi2) (15) 

A number of later surveys confirm this conclusion (§4, Peebles 1979, DP83, Bean et al 
1983, Hale-Sutton et al. 1989, Fisher et al. 1994b). 

Velocity distributions derived from N-body simulations vary. Efstathiou et al. (1988) 
note that the distribution of pairwise velocities in their simulations differs substantially 
from an exponential. The difference is particularly noticeable at small relative velocities, 
where P{vi2) is much fiatter than the exponential distribution. Fisher et al. (1994b) 
confirm this behavior in the simulations they use to interpret their data. 

Following both coUisionless and baryonic particles in their N-body simulations. Gen 
and Ostriker (1993) find that the distribution of 3-d szra^f/e- "galaxy" peculiar velocities is 
well represented by an isotropic exponential: 

/(vi)dV =exp-/9|vi|d3vi (16) 

where /? = \l2j<T\. For this distribution, the probability of measuring a velocity with 
magnitude v\ is proportional to v\e^~^^^^y 
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We would like to explore not only the moments of P{vi2) but also the functional form 
of the distribution. We can derive a fitting function for P{vi2) to compare with the N- 
body model in Cen and Ostriker (1993) as follows. If the velocity field is isotropic, the 
probability of measuring one component v" of a single-galaxy velocity is 

/oo /"OO ; 
/ exp [-^\^{vfr + {vlr + {vfr\] dvfdvf 
-co J — CO 

= (l<l + ;^)exp[-/9|<|] (17) 

The second moment of this one-dimensional, szra^f/e-galaxy distribution is a/2(Ji. A similar 
transformation yields the distribution of one component of the pairwise velocity P{vi2), 
where — ^"(1) ~ ^?(2)- To simplify the notation, we define the variables Ui — v"{i) so 
that — ui — U2. Assuming the velocities of the galaxies are isotropic and independent, 

/oo 
Fl«2 + U2)Fl(u2)du2 
-00 

yoo II 
= J_JM2 + + ^)(h2| + ^)exp [-[^{^2 + + \U2\)] du2 

= + ^ + ^(l<2l + ^)\ exp (-/9|<2l) (18) 

The second moment of this one-dimensional pairwise distribution P{vi2) is 2(Ji, i.e. twice 
the 3-d single-galaxy velocity dispersion. 

Figure 1 shows the normalized distributions P{vi2) given by Equations 15 and 18. The 
distributions are normalized to have the same second moment; (J12 = 500kms~^ for the 
distribution in Equation 15 and (Ji = 250kms~^ for the distribution in Equation 18. If 
P{vi) is exponential and we mistakenly assume that P{vi2) is exponential with (J12 = 2(Ji, 
we overestimate the number of pairs at small separations. It is important to emphasize 
that an exponential distribution of szra^f/e-galaxy velocities does not imply an exponential 
distribution of pairwise velocities. 

The distribution P{vi2) in Equation 18 is qualitatively similar to pairwise velocity 
distributions measured directly from N-body simulations. Efstathiou et al. (1988) and 
Fisher et al. (1994b) plot these distributions for a range of initial conditions and epochs, 
and in each case, the core of the distribution is noticeably fiatter than exponential. This 
small-scale fiattening may be partially responsible for the discrepancy observed by Fisher 
et al. (1994b) between their isotropic exponential model for C(rp,7r) and the CF calcu- 
lated directly from their N-body simulations. Of course, as Fisher et al. (1994b) and Cen 
and Ostriker (1993) note, pairwise velocities in the real universe and in simulations are 
anisotropic and correlated; thus the assumptions that go into the derivation of Equation 
18 are not entirely valid. Although Equation 18 only approximates the behavior of V12 in 
the simulations, it provides a useful guide to departures from the exponential distribution 
in Equation 15. 

To summarize, we test the underlying velocity field by fitting models to C(rp, vr) using 
both exponential pairwise distributions and exponential szra^f/e-galaxy distributions. For 
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the latter case, we use the pairwise distribution in Equation 18 for the convolution Equation 
8. We test the validity of Equation 18 in §4.2. We emphasize that neither distribution is 
based on a physical model; Equation 15 is motivated by observations, and Equation 18 is 
a simple empirical model which ignores correlations and anisotropy in the velocity field. 
Clearly, a physical model for the velocity distribution function in the non-linear regime 
would be extremely useful (e.g. Kofman 1991 and references therein). 



3.6 Scale Dependence of the Velocity Moments 

If clustering is statistically stable, the cosmic virial theorem implies that a 12 scales as 
r^~'^ where 7 is the slope of the spatial CF. Thus for 7 in the observed range (~ 1.8-1.9), 
(712 is nearly constant with scale. On large scales, linear theory requires that (J12 decrease 
slowly with scale (Peebles 1980). The scaling of (J12 in N-body simulations varies according 
to the details of the model and the algorithm. Fisher et al. (1994b) adjust the scaling of (J12 
from their N-body simulations to match the observed CF and to reproduce both the cosmic 
virial theorem on small scales and the linear theory prediction on large scales. Because 
the behavior of (J12 derived by Fisher et al. is essentially scale-independent, and because 
we don't know what resemblance the simulations bear to the real universe, we choose to 
ignore the scale dependence of (J12. We show in §4 that this choice does not hinder our 
ability to measure a 12 on small scales even in simulations where a 12 varies significantly 
with scale. 

Our knowledge of the scaling of the first moment, ^12, is similarly vague. Because (J12 is 
likely to be significantly larger than V12 on small scales, we choose to calculate (J12 without 
modeling the infall pattern. To investigate the effects of infall on our derivation of (J12, we 
adopt the model used by DP83 derived from the scale-invariant solution to the truncated 
BBGKY heirarchy (Davis k Peebles 1977): 

vi2{r) = f^T^ (19) 

^ ^ l + (r/ro)2 ^ ^ 

The similarity solution implies F — 1 (DP83). 



3.7 Model Fitting 

Standard analysis is inadequate for correlation analysis for two reasons. First, 
estimates of the CF, ^, at different separations are correlated. The correlation between 
points effectively reduces the number of degrees of freedom; if correlations are large, they 
can significantly alter the fit. Second, Fisher et al. (1994a) report that the distribution of 
^ is significantly non-normal over portions of the (rp, 7r)-plane. 

In order to address the covariance between ^ at different separations, we use techniques 
of principal component analysis elucidated by Fisher et al. (1994a). We compute the 
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covariance matrix using 100 bootstrap resamplings of the original dataset. Because the 
covariance matrix of ^ is real and symmetric, it can be diagonalized by an orthogonal 
transformation. This transformation defines a set of linear combinations of the original 
data points which are linearly independent and are therefore appropriate for analysis. 
The transformed statistic is 

~2 ~ ii, model) (0C\\ 

X — ^2 V^^/ 

i-l 

where ai are the eigenvalues of the original covariance matrix and and model '^^^ 
the transformed data and model vectors, respectively. The columns of the matrix which 
transforms ^ to ^ are simply the eigenvectors of the original covariance matrix V. 

The non-normality of ^ is more problematic. Fisher et al (1994a) use N-body sim- 
ulations to show that this distribution is slightly but significantly skew on the scales of 
interest. Thus the standard distribution does not accurately represent the expected 
distribution of our statistic. Although this non-Gaussianity mitigates some of the sta- 
tistical advantages of principal component analysis, we use to address the covariance 
between points and quote the confidence interval under the assumption of Gaussian er- 
rors. Because the error distribution is likely to be broader than a Gaussian, we interpret 
borderline rejections ~ .01) cautiously. 
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4 TESTS OF THE METHOD 



In this section, we test the accuracy of the procedure outUned in §3 using simulated 
data. M. Crone and A. Evrard kindly provided an N-body simulation for this analysis. 
Crone and Evrard (1994) follow the motion of 262,144 CDM particles in a 64 h'^Mpc 
box using a particle-particle/particle-mesh algorithm. Initial density perturbations have a 
power-law spectrum 6(k.) — \k\"' where n — —1, and the cosmology is flat {il — 1,A = 0). 
Using a random subsample of the particles, we construct mock redshift surveys by viewing 
the particles from corners of the box and assigning a uniform Hubble flow. From these 
mock surveys, we compute C(rp, vr), w{rp), and ^(r) as we would for the real data. Because 
the simulation volume is small, we do not sample enough independent structures to draw 
large numbers of mock surveys; this restriction limits our ability to test for bias in the 
techniques. Still, the simulation serves as a convenient and physically motivated set of 
known phase-space coordinates. We do not intend to test these models against the data; 
we use them only to establish the accuracy and the limitations of our procedure. 



4.1 ^(rpjir) from the Simulation 

The velocity held of the N-body simulation is very hot; (Ji2(0.5/i~^ Mpc) = 1435 kms~^. 
Because we cannot measure C(rp,7r) accurately beyond 30 h~^M.]yz, we cannot test our 
procedure adequately with the raw velocity held. Following Davis et al (1982), we cool 
the velocity held by dividing all particle velocities by a factor of three. This procedure is 
roughly (though not precisely) equivalent to lowering the mean mass density to il ~ 0.11 
(Davis et al 1982). The reduced dispersion is ai2{0.bh~^ Mpc) ~ 500 km s"-*-. 

Figure 2a shows C(rp,7r) calculated from a mock survey of the simulated galaxy dis- 
tribution. The elongation along the line of sight caused by peculiar velocities is clearly 
evident. Compression of the contours at large Vp hints at the infall pattern discussed in 
§3.2. Figure 2b shows C(rp,7r) in the absence of peculiar velocities. Although the spatial 
distribution is nearly isotropic on scales less than the correlation length (~ 5/i~^Mpc), 
the shape of the structures themselves distort C(rp, vr) on larger scales. Because we ascribe 
anisotropy in C(rp, vr) to peculiar velocities, this intrinsic spatial anisotropy fundamentally 
limits our ability to extract the velocity moments (particularly ^12) from flnite redshift 
surveys (Fisher et al 1994b, Marzke et al 1994b). 



4.2 Inversion of w{rp) for ^{r) 

Figure 3 shows how the observational cutoff iimax affects the deprojection of ^(r). The 
solid line is a pure power law. We use this power law to make an analytic model of ^(rp, vr) 
using Equation 8 and then calculate w{rp) numerically. We then invert w{rp) for ^(r) 
using Equation 12. We show inversions for i^max — 30/i~"'" Mpc and for iTmax — QOh~^ Mpc. 
Deviations from a power law are signiflcant beyond 10 h~^Mpc for iVmax — 30/i~^Mpc. 
We therefore restrict our analysis of ^(r) to r < 10h~^ Mpc. 
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Figure 4 shows the fit to w{rp) measured from the simulation. Figure 5 compares ^(r) 
computed from the mock survey with the CF calculated directly from pair counts using 
the known spatial coordinates. To the extent that the simulated CF can be represented 
by a power law, the inversion of w{rp) reproduces the spatial CF on scales smaller than 
~10 Mpc. 



4.3 Recovering a Known Velocity Field 

To test our ability to recover a known velocity field, we start with a simple model. 
For each particle in the N-body simulation, we replace the computed velocities with 
random deviates drawn from the isotropic exponential distribution in Equation 16 with 
(Ji = 250kms~^. We preserve the clustered spatial distribution; we replace only the 
computed velocities. Figure 6 shows the distribution of 1-d pairwise velocities measured 
from the simulation along with the prediction of Equation 18. Equation 18 predicts the 
measured pairwise distributions well. 

Figure 2c shows C(rp,7r) calculated from a mock redshift survey of this model. On 
large scales, the characteristic compression caused by infall, which is subtle but clear in 
figure 2a, is clearly absent in figure 2(c) where V12 — 0. Figure 7 shows fits to C(rp,7r) 
for the first six broad bins in Vp. The solid line represents the convolution in Equation 8 
with P{vi2) given by Equation 18. The fit is good. We recover the input ai to within one 
standard deviation: cri{rp < O.Sh~^ Mpc) = 260 it Hkms"-*-. 



4.4 Recovering (J12 from the N-body Simulation 

In the N-body simulation, the velocity moments vary with pair separation (cf. Efs- 
tathiou et al. 1988, Fisher et al. 1994a). Figure 8 shows (J12 and V12 as a function of pair 
separation. We measure these moments directly from the particle positions and velocities 
computed in the simulation. Because we assume in our model that the velocity dispersion 
is independent of scale, we worried that the dispersion (J12 we infer from ^{rp, it) might be 
biased. The open squares in Figure 8 show (J12 measured from C(rp,7r) using the model 
(Equation 8) and ignoring infall. The measured dispersion follows the true dispersion 
remarkably well. With only one exception, the measured a 12 falls within one standard 
deviation of the true dispersion. 

The dotted line in Figure 8 shows the measured a 12 assuming the infall model in 
Equation 19 with F — l.b (this model is shown as the dashed line in Figure 8). The 
inclusion of infall actually degrades the measurement of (J12. Although the dispersion 
at Tp < l/i~^Mpc is consistent with the true dispersion, the measurements at larger Tp 
are biased high. DP83 note a similar discrepancy in their analysis of the simulations 
described by Efstathiou and Eastwood (1981). The source of this discrepancy is unclear, 
but there are several possibile culprits: the spatial CF deviates from a power law in the 
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N-body simulation, and in detail, (J12 is neither isotropic nor scale-independent. Because 
we recover (J12 within one standard deviation without including infall, we derive (J12 for 
the data using the no-infall model. We show results for (J12 including infall to emphasize 
that the dependence of (J12 on Vp may reflect the scaling of V12 more than it reflects the 
true behavior of (J12. 
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5 RESULTS 



5.1 ^(rpjir) for the Data 

Figures 9 and 10 show contours of ^(rp, vr) for CfA2N, CfA2S, SSRS2 and CfA2+SSRS2. 
In order to facilitate comparisons with other surveys and to depict more clearly the dis- 
tortions in C(rp, it), we display ^{rp, it) with linear lh~^ Mpc binning. We emphasize that 
this binning is used for display purposes only; for model fitting, we use logarithmic binning 
exclusively. 

Pairs are MV weighted in Figure 9 and uniformly weighted in Figure 10. The choice 
of weighting significantly affects C(rp,7r). Although the sensitivity to weighting is small 
in CfA2+SSRS2, the remaining fiuctuations emphasize that surveys on this scale do not 
probe a "fair" volume of the universe. The fiuctuations in C(rp,7r) with weighting are not 
surprising given the redshift distributions of these samples (Geller and Huchra 1989, Park 
et al 1994, da Costa et al 1994). In CfA2N, uniform weighting accentuates the generally 
low-density features nearby. Uniform weighting also increases the contribution of faint 
galaxies, which are less clustered than brighter galaxies (Giovanelli & Haynes 1989,Vogeley 
et al 1994). On the other hand, MV weighting enhances the contribution of the Great 
Wall, which contains most of the groups and clusters in CfA2N. The clustering amplitude is 
consequently larger for MV weighting , and the orientation of this huge structure elongates 
^(rp,7r) in the Tp direction. 

CfA2S reveals a similar trend. In this sample, the Perseus-Pisces wall at ~ 5, 000 kms~^ 
lies near the peak of the distribution 4'{r)V{r) where (f) is the selection function and V is 
the volume of the shell at distance r. In this case, MV weighting enhances the low-density 
regions in the foreground and background, lowering the amplitude of ^(rp,7r). The shape 
of ^(rp, tt) on small scales is only mildly sensitive to the weighting in CfA2S. Like the Great 
Wall in CfA2N, the Perseus-Pisces wall runs perpendicular to the line of sight; as expected 
from the redshift distribution, uniform weighting enhances correlations in the Tp direction. 

The SSRS2 is fairly insensitive to the weighting scheme on scales smaller than ~ 10h~^ Mpc. 
This insensitivity is reasonable, as the most obvious structure crosses the survey at an angle 
to both the Tp and tt axes. On larger scales, MV weighting strongly enhances line-of-sight 
correlations, again refiecting the small number of independent structures sampled in this 
volume. Although uniform weighting provides insight into the fiuctuations in C(rp,7r) 
within the sample, it is a somewhat unnatural and inefficient weighting scheme; we there- 
fore drop the uniformly weighted CF from further consideration. We calculate the spatial 
CF and the pairwise dispersion using MV weighting exclusively. 



5.2 Spatial Correlation Function 

Figure 11 shows power-law fits to w{rp) for each sample. We summarize these fits 
in Table 3, which lists the parameters of the spatial CF obtained from the inversion in 
Equation 12. As discussed in §3, we transform the data vector w{rp) to diagonalize the 
covariance matrix and calculate from Equation 20. A single power law fits w{rp) for 
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each sample at separations smaller than 10 h Mpc. The correlation lengths are roughly 
consistent for CfA2N and SSRS2: tq = 5.83 ± 0.18 and 5.08 ± 0.23, respectively. The 
slopes are also consistent: for CfA2N, j = -1.80 ± 0.03; for the SSRS2, j = -1.85 ± 0.06. 
In CfA2S, the amplitude of the clustering is significantly weaker and the slope is steeper: 
To = 4.75 lb 0.20h~^ Mpc and 7 = —1.99 it 0.05. Because tq and 7 are strongly correlated, 
these differences are not as large as they appear; if we constrain 7 to be -1.8, tq = 5.4 for 
CfA2S. For consistency, we use the fitted values for both the amplitude and the slope to 
calculate (J12. 



5.3 Pairwise Dispersion for CfA2 and SSRS2 

Figures 12-15 show cuts of ^(rp, vr) at each rp for CfA2N, CfA2S, SSRS2 and CfA2+SSRS2 
respectively. We include fits of Equation 8 based on the two distributions P{vi2) in Equa- 
tions 15 and 18. Dashed lines represent the single-galaxy exponential (Equation 18); solid 
lines represent the pairwise exponential distribution. Again, we neglect the possible con- 
tribution of the mean pairwise infall V12. 

Error bars in Figures 12-15 represent the raw bootstrap variances. Tables 4-7 list 
estimates of (J12 at each value of rp for each sample. Columns 3-5 show fits to the single- 
galaxy exponential in Equation 18. Columns 6-8 give the fits for a pairwise exponential. 
Tables 4-7 also include fits using the infall model in Equation 19 with F=l. 

In order to facilitate comparisons with models and with other surveys, we compute 
a best estimate for (J12 as the mean of the first three bins in rp. This mean, cri2{rp < 
0.8/i~^ Mpc), is listed for each sample in Table 8 along with its associated error. The 
variation in cri2(< 0.8) among samples is striking. For the exponential P{vi2) with no 
infall, (7i2(< 0.8) = 647±52kms-^ 367±38kms-^ and 272±42kms-^ for CfA2N, CfA2S 
and SSRS2, respectively. In general, including infall increases the estimate of cri2(< 0.8), 
consistent with DP83. 

The behavior of (J12 as a function of rp varies with the infall model. Figure 16 shows 
Ci2(?'p) for each sample. If we include infall, we find that (J12 generally rises slowly or 
is effectively constant with rp, consistent with the observations of DP83 and Mo et al. 
(1993). However, if we ignore infall, the dispersion drops with rp in each case. Because 
of the difficulties with the infall model at large rp highlighted in §4, we suggest that the 
dependence of (J12 on rp measured with this technique should be interpreted with caution. 

Using the individual samples with the standard binning, the errors preclude a strong 
rejection of either model for the velocity distribution function. In order to increase the 
signal-to-noise at small rp, we recompute C(rp,7r) for the combined CfA2+SSRS2 using 
linear, lh~^ Mpc bins (see figure 9). Figure 17 shows the cut in C(rp, vr) at rp < lh~^ Mpc 
along with fits based on the velocity distributions in Equations 15 and 18. The single- 
galaxy exponential, represented by the dashed line, is clearly incompatible with the data 
(x^/d.o.f. = 37.7/13,(5 = 2 X 10"'*). The single-galaxy exponential is too fiat at small 
velocities and drops too quickly at large velocities. On the other hand, the pairwise expo- 
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nential reproduces the data well (x^/d.o.f. = 10.5/13,(5 = 0.65), and the corresponding 
dispersion 565ib33 is consistent with our best estimate derived above, cri2(< 0.8) = 540ib36. 
Thus the velocity field in the nearest ~ 120/i~^Mpc is remarkably consistent with an 
isotropic exponential distribution of pa«™«se velocities, consistent with previous surveys ( 
e.g., DP83, Mo et al 1993, Fisher et al 1994b). 

Table 8 summarizes the most important results of this analysis: the pairwise disper- 
sions cri2(< 0.8) measured from the CfA2 and SSRS2 samples. For the combined sample, 
we estimate the error in (J12 from the scatter in the individual sample estimates. This er- 
ror (labeled "external") is simply the weighted sample standard deviation. Two points are 
clear: 1) (J12 varies broadly among volumes < 10'^ Mpc^, and 2) (J12 for CfA2+SSRS2 
significantly exceeds previous estimates of the pairwise velocity dispersion, which yield 
(712 ~ 300. The agreement on (J12 between Fisher et al. (1994b) and DP83 may be for- 
tuitous; CfAl fully samples only one rich cluster, and the IRAS survey is dominated by 
spiral galaxies, which avoid the cores of clusters. Indeed, Mo et al. (1993) and Zurek et al. 
(1994) find that the value of (J12 derived from CfAl varies significantly depending on how 
the velocity field around the Virgo cluster is treated. With no corrections to the velocity 
field, the estimate of (J12 for CfAl derived by Zurek et al. (1994) is consistent with (J12 
derived from CfA2+SSRS2. Furthermore, (J12 for ellipticals and SO's is significantly larger 
than it is for spirals (Mo et al. 1993, Marzke et al. 1994b). Thus the variations in the 
observed (J12 appear to stem from the contributions of a few rich clusters. We test this 
hypothesis in the following sections. 



5.4 The Cluster Contribution 

To guage the contribution to a 12 from rich clusters, we mask out galaxies within three 
degrees of Abell clusters with richness class R > 1. Figure 18 depicts C(rp,7r) for these 
pruned samples. Not surprisingly, the removal of clusters suppresses the extension of 
^(rp,7r) along the vr direction. 

Table 3 includes fits to w{rp) with clusters removed. Except in CfA2S, the slope of 
the spatial CF does not change when we remove clusters, but the amplitude decreases 
significantly. In CfA2S, removing clusters increases the amplitude and fiattens the slope. 
These changes in ^(r) have little effect on the calculation of (J12 on small scales. 

Tables 4-7 list a 12 derived from the pruned samples. In most cases, a 12 drops signifi- 
cantly when we remove clusters. The change is most dramatic in CfA2N, where cri2(< 0.8) 
drops from 647 it 52 to 313 it 39. Excluding all galaxies within 2500 kms~^ in CfA2N does 
not change this result significantly; in this case, (J12 = 327 it 41. Thus the infall pattern 
around the Virgo cluster does not affect our derivation of a 12. The cores of rich clusters 
are the dominant source of variation in a 12. 

In CfA2S, we measure a similar decline in (J12 from 367 it 38 km s~^ to 264 it 38 km s~^ 
when we remove clusters. On the other hand, the pruned SSRS2 sample is consistent 
with the SSRS2 as a whole. Because the SSRS2 contains no very rich clusters like Coma, 



19 



this insensitivity is not surprising (Ramella et al 1995). For the combined sample, we 
find cri2(< 0.8) = 295 it Slkms"-*- after removing clusters. Thus the difference between 
CfA2+SSRS2 and DP83 can be attributed entirely to the presence of rich clusters in CfA2. 
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6 DISCUSSION 



6.1 Comparison with Cosmological Models 

The pairwise velocity dispersion measured by DP83 has stood for a decade as one of 
the strongest arguments against the standard CDM paradigm. Taken at face value, our 
results appear to alleviate the conflict. Figure 19 displays (J12 measured from CfA2+SSRS2 
along with a compendium of model predictions. We include only variants of CDM; we 
do not consider possible mixtures of CDM with other species, which provide an equally 
good match to the measured power spectrum (Summers et al 1994, Klypin et al 1993). 
Although unbiased flat CDM still fails to match the data, biased il — 1 models are no 
longer ruled out by the pairwise velocity dispersion. 

Our data also appear consistent with the open models discussed by Davis et al (1985). 
Interestingly, open models seem to provide the best match to a wealth of observational 
constraints: the measured power spectrum, void probability function, and small-scale ve- 
locity held are all consistent with Qh 0.2 CDM (Vogeley et al. 1992, 1994, Park et al. 
1994, Fisher et al. 1992). On the other hand, biased CDM, which satisfles our constraint 
from (J12, does not reproduce the observed void probability function (Vogeley et al. 1995). 

Unfortunately, the interpretation of (J12 is complicated by the presence of rich clusters. 
Because we and others have shown that (J12 is very sensitive to the exclusion of clusters, the 
pairwise dispersion we measure cannot be construed as the true "held" velocity dispersion. 
The open squares in Figure 18 show (J12 from CfA2+SSRS2 after removing Abell clusters 
with R > 1. If the model predictions represent the field, then the constraints on all variants 
of CDM still stand. 

Clearly, the deflnition of the "held" poses a signiflcant problem for the interpretation 
of the small-scale velocities in both the models and in the real universe. For example, 
Couchman and Carlberg (1992) find that (J12 varies between 100 kms"-*- and 1000 kms"-*- 
even in moderately dense regions. The value for (J12 quoted in Couchman and Carlberg 
(1992) and reproduced in Figure 19 does not include the richest cluster in their simulation. 
With the current limitations on both observations and theory, the effort to constrain 
cosmological models demands careful, consistent flltering of both simulations and data to 
establish fair grounds for comparison. 

In any clustering scheme, gravity builds up a continuum of bound systems. The removal 
of clusters is somewhat arbitrary, particularly when the surrounding large-scale structure 
is complex and the observed (or simulated) volumes are small. Ideally, one would like 
the simulations to reproduce (J12 averaged over all galaxies. On the theoretical side, this 
demand calls for a detailed treatment of the interaction among galaxies in a large enough 
volume to include an ensemble of large, dense structures. Observationally, the requirement 
is a fair sample of the universe. 

These are tall orders. In order to size up the task before us, we ask the following 
questions : 1) how does the distribution of cluster velocity dispersions affect (J12 and 2) on 
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what scale does (J12 converge to its global mean? 



6.2 Pairwise Dispersion from Virialized Systems 

We develop a simple model to describe the variation of (J12 with survey volume. Red- 
shift surveys indicate that a substantial fraction of galaxies aggregate into bound sys- 
tems with velocity dispersions ranging from ~ lOOkms"^ in loose groups to more than 
1000 km s~^ in rich clusters (Geller and Huchra 1983, Nolthenius and White 1987, Ramella 
et al 1989, Moore et al 1993, Zabludoff et al 1993). Because we know from §5.4 that 
dense systems contribute significantly to the field velocity dispersion (J12, we can predict 
a 12 and its variance in a given volume from the distribution of group and cluster velocity 
dispersions, n{a) (cf. Carlberg 1994). We can obtain n{a) either from the measured dis- 
tributions (Zabludoff et al 1993, Frenk et al 1991) or directly from the power spectrum of 
fiuctuations using the Press-Schechter formulation for the evolution of the mass spectrum 
(Press and Schechter 1974). 

By definition, the statistic (J12 is pair- weighted. Because the number of pairs in a 
cluster is proportional to the square of the number density of galaxies within the cluster, /o, 
the mean pairwise dispersion is 



where (Jmrn and (Jmair limit the range of group and cluster velocity dispersions under con- 
sideration. The largest measured dispersion is cr = 1436 km s"-*- (Oegerle 1994, private 
communication); extrapolation of n{cr) to larger values is uncertain. We assume in Equa- 
tion 21 that both members of a pair inhabit the same system, a reasonable assumption on 
scales < l/i~^Mpc. Dropping this assumption increases the expected (J12. The relation 
between the number density of galaxies within a single cluster and the cluster velocity 
dispersion, p{cr), requires an assumption about the distribution function. We parameterize 
this dependence as /o oc cr™': 



) ^ mm ^ 

r'Jr^a^ (j2f»n((j) da 

mm ^ 



z _ 2iii2^i£ (22) 

12 ^ rl^max _2m„/_\ J_ \^^/ 



For an isothermal sphere, m — 2. Because observational constraints on m are weak, we 
assume m — 2 throughout this analysis. If m > 2, both (J12 and the predicted variance in 
(712 increase. We ignore the contribution of an unbound, true field population. Including a 
low- dispersion field would reduce the estimated (J12 but would not contribute significantly 
to the variance in a 12. 
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6.2.1 (712 from the Observed Distribution of Velocity Dispersions 

ZabludofF et al. (1993) measure the distribution of velocity dispersions for dense sys- 
tems. Their sample includes dense Abell clusters with richness class > 1 as well as 
a sample of dense groups selected from the first four slices of CfA2(Geller & Huchra 
1989,Huchra et al. 1994). The authors conclude that the number density of systems de- 
clines exponentially with velocity dispersion: n{a) oc iQ-^-^^^^<^ _ Using a similar sample 
of groups from CfAl (Huchra et al. 1983), Moore et al. 1993 find a steeper slope for n{a), 
roughly proportional to iQ-^-^^'^<^ _ Their comparison with the distribution of cluster ve- 
locity dispersions culled from the literature (Frenk et al. 1990) suggests a discontinuity 
in n{a) at cr ~ 700kms~^, where the distribution becomes fiatter. The authors attribute 
this fiattening to field contamination in the cluster redshift distributions, which artificially 
infiate the measured dispersions. The n{a) inferred from the distribution of cluster X-ray 
temperatures falls somewhere in between the results of Zabludoff et al. (1993) and Moore 
et al. (1993); The X-ray temperature distributions follow a power law corresponding to 
n{a) oc over the observed range (Edge et al. 1990, Henry & Arnaud 1991). 

We use these measurements of n{a) to estimate the volume-averaged (J12 according to 
Equation 22. We show the range of n{a) from the observations as dashed lines in Figure 
20. Following convention, we plot the cumulative distribution n{> a). Each distribution 
is normalized to the mean density of Abell clusters with R > 1 n{a > 700kms~^) = 
6.6 X lO-^h^Mpc-^ (Zabludofl: et al. 1993). 

Figure 21 shows the expected mean pairwise dispersion,(Ji2, as a function of Umax from 
Equation 22 for the observed range of n{a). The convergence of (J12 depends on the shape 
of n{a) and on the extrapolation of the observed distributions beyond the observed range 
of a. For steeper n{a), (J12 converges to lower values and at lower Umax- If we simply 
extrapolate the distribution n{a) derived by Zabludoff et al. (1993), ai2 does not converge 
until 

'^max ~ 3500 km s well beyond the observed range. 

The velocity dispersion Uconv at which (J12 converges determines a minimum volume 
required for the statistic to be representative. In order to predict the variance of (J12 
for a given survey volume V , we perform a simple Monte-Carlo simulation. The number 
of clusters in the hypothetical sample volume V is U[n(> Umin) — ri{> Umax}]- As 
in Zabludoff et al. (1993), we assume Umin — 100 km s"-*-. The cutoff Umax is a free 
parameter. At each volume, we construct several Monte-Carlo realizations of a sample of 
clusters drawn from the distribution n(cr) . We then measure (J12 = YliLi o'f"^'^'^ / YliLi '^1"^ ■ 
Because of the large weight given to high-dispersion systems, our analysis of the bias and 
variance in (J12 is much more sensitive to (Jmair than it is to Umin- 

Figure 22 shows the distribution of the fractional error in cri2,^ci2/c''i2, over the range 
V — 10^ — 10^/i~^Mpc^ for the range of n{a) derived from observations. Arrows indicate 
the volumes of CfAl, CfA2N and CfA2+SSRS2. The solid line represents the median (J12; 
dotted lines enclose 68% of the results for each volume. Because (J12 is weighted by the 
square of the density, the shot noise caused by the one or two richest clusters in each sample 
dominates the variance at small volumes. The distribution of a 12 is noticeably skew with a 
tail toward large a 12. Furthermore, a 12 is biased toward low values in samples with small 
volumes. As the volume increases, samples include significant numbers of clusters with 
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c ~ o'conv, and the variance drops accordingly. On the scales of the surveys, the variance 
and bias are larger than the random errors in (J12, especially for the shallow distribution 
in Zabludoff et al (1993). 

Figure 22 also shows the sensitivity to (Jmair- Increasing (Jmair from 2000 kms"-*- to 
3000 kms"-*- affects the convergence of (J12 much less than changing the slope of n{cr) over 
the observed range. As expected, extrapolating to larger a slows the convergence of (J12 
with volume. There is no reason to expect that the observed distributions n{cr) extend 
beyond the observed range; in the next section, we investigate a more physical model for 
the extremes of n{cr). 



6.2.2 (712 from the Press-Schechter n{a) 

The mass distribution of virialized systems at any epoch is largely determined by the 
low-order moments of the initial density perturbations. Assuming spherical collapse of 
Gaussian perturbations. Press and Schechter (1974) describe the evolution of the distribu- 
tion of mass clumps for various assumptions about the power spectrum of initial density 
perturbations. Although the ingredients of the model are approximate, N-body simulations 
of the non-linear (and non-spherical) collapse of structure show a remarkable agreement 
between the Press-Schechter prediction and the measured n{M, z) for a wide range of ini- 
tial conditions (Efstathiou et al. 1988, Carlberg and Couchman 1989). Recent attempts 
to improve the analytic model suggest that the low and high-mass ends of n{M, z) may 
not be accurate (e.g. Peacock and Heavens 1990, Monaco 1994), but we are primarily 
interested in the regions in between. 

To estimate the distribution of cluster velocity dispersions, we calculate the Press- 
Schechter distribution in the form presented by Carlberg et al. (1994): 

47rV27r (J* rflncr ■ ' 

where v — 1.68(l + z)/(cr8 A(M)), erg is the normalization of the mass fluctuation spectrum, 
and A(M(r)) is the rms mass fluctuation on the scale r. We use the analytic fit to A(r) 
for fiat CDM given by Narayan and White (1988). The constant defines the relation 
between a and mass M: (7 = CyHoR{l + zfl"^, where M = 4.%poR^/'i. White and Frenk 
(1991) measure — 1.18 for the spherical collapse of a tophat sphere. For standard CDM, 
the only free parameter is the normalization, erg. 

Solid lines in Figure 20 show the cumulative distribution of velocity dispersions given 
by the Press-Schechter prescription for two normalizations of the power spectrum. The 
normalization cr^ — 1.0 corresponds approximately to the measured quadrupole anisotropy 
of the microwave background (Wright et al. 1992). The second normalization, erg = 0.5, 
indicates the low end of the range currently favored by models of biased galaxy formation 
(Couchman and Carlberg 1992, Cen and Ostriker 1993). The distribution of Abell cluster 
velocity dispersions (Zabludoff et al. 1993) is significantly shallower than the CDM predic- 
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tions for any acceptable value of erg (this discrepancy is unlikely to be caused entirely by 
projection effects; c.f. Zabludoff et al 1990 and Frenk et al 1990). 

Figure 21 shows the convergence of (J12 computed from the models. Even for erg = 1.0, 
(J12 converges within the observed range of cluster velocity dispersions. (J12 converges 
to 1239 kms~^ for COBE-normalized CDM, remarkably consistent with the dark-matter 
velocity dispersion derived by Cen and Ostriker (1993) for ~ lh~^ Mpc pair separations. 

Figure 23 shows 5(Ti2/(Ti2 given by the Monte-Carlo procedure described above but 
using n{a) from the models. As expected, the convergence of (J12 with volume depends 
strongly on erg. For erg = 0.5, a 10% measurement of (J12 is obtained in a survey volume 
~ lO^Mpc^. For unbiased CDM, the required survey volume is much larger, approximately 
5 X lO'^Mpc^, which is larger than any existing redshift survey but comparable with the 
planned Sloan Digital Sky Survey (Gunn and Weinberg 1995). 

We conclude that the variation in (J12 among subsamples of CfA2+SSRS2 is consis- 
tent with expectations based on the observed distribution of group and cluster velocity 
dispersions. The volume required for a 10% measurement of (J12 is larger than the volume 
surveyed by CfA2+SSRS2 in an unbiased, flat CDM universe. Furthermore, we show that 
(J12 is biased low by approximately 20% on the scale of CfAl (used by DP83) for the 
same model. The bias and variance are worse if n{a) observed by Zabludoff et al (1993) 
extends to larger a. Convergence occurs on a smaller scale for high-bias CDM models. 
Our model is admittedly approximate; cluster-cluster correlations increase the expected 
variance, while velocity bias tends to decrease it. It is clear, however, that large redshift 
surveys as well as large, high-resolution simulations are required for a robust comparison 
between observed and predicted pairwise velocity dispersions. 
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7 CONCLUSIONS 



We combine the CfA Redshift Survey (CfA2) and the Southern Sky Redshift Survey 
(SSRS2) to estimate the pairwise velocity dispersion of galaxies on a scale of ~ lh~^ Mpc. 
Both surveys are complete to an apparent magnitude limit B{0) — 15.5. Our sample 
includes 12,812 galaxies distributed in a volume 1.8 x 10'^ Mpc^. We conclude: 

1) The pairwise velocity dispersion of galaxies in the combined CfA2+SSRS2 redshift 
survey is (J12 = 540 km s"-*- it 180 km s"-*-. Both the estimate and the error in (J12 exceed 
the canonical values (J12 = 340 it 40 measured by Davis and Peebles (1983) using CfAl. 

2) We derive the uncertainty in (J12 from the variation among subsamples with volumes 
on the order of 7 x 10^ Mpc^. This variation is nearly an order of magnitude larger than 
the formal error 36 kms"-*- derived from least-squares fits to the CfA2+SSRS2 correlation 
function. This variation among samples is consistent with the conclusions of Mo et al 

(1993) for a number of smaller surveys and with the analysis of CfAl by Zurek et al 

(1994) . 

3) When we remove Abell clusters with R > 1 from our sample, the pairwise velocity 
dispersion of the remaining galaxies drops to 295kms~-'-ib99kms~-'-. This dispersion agrees 
with the results for CfAl (DP83) and with the results for the IRAS 1.2 Jy survey (Fisher 
et al 1994b), both of which are biased against rich clusters. 

4) The dominant source of variance in a 12 is the shot noise contributed by dense 
virialized systems. Because (J12 is pair- weighted, the statistic is sensitive to the few richest 
systems in the volume. This sensitivity has two consequences. First, a 12 is biased low 
in small volumes, where the number of clusters is small. Second, we can estimate the 
variance in a 12 as a function of survey volume from the distribution of cluster and group 
velocity dispersions n{a). For either a COBE-normalized CDM universe or for the observed 
distribution of Abell cluster velocity dispersions, the volume required for (J12 to converge 
(^0-12/0-12 < 0.1) is ~ 5 X lO^/i-3 Mpc^, larger than the volume of CfA2+SSRS2. 

5) The distribution of pairwise velocities is consistent with an isotropic exponential 
with velocity dispersion independent of scale. The inferred szra^f/e-galaxy velocity distribu- 
tion function is incompatible with an isotropic exponential. Thus the observed kinematics 
of galaxies differ from those measured in the hydrodynamic simulations of Cen and Ostriker 
(1993). On the other hand, our observations appear to be consistent with the distribu- 
tion function measured by Zurek et al. (1994) using coUisionless simulations with much 
higher resolution than Cen & Ostriker (1993). The large dynamic range in the simulations 
by Zurek et al. (1994) affords an accurate treatment of galaxy interactions, which dra- 
matically alter the dynamical evolution the galaxy distribution (Couchman and Carlberg 
1992, Zurek et al. 1994). The agreement between the observed velocity distribution and 
the one predicted by these high-resolution simulations may be another clue that mergers 
play an important role in the evolution of galaxies and large-scale structure. 

We thank Mary Crone and Gus Evrard for generously providing one of their N-body 
simulations for the tests in §4. We also thank the referee, Jim Peebles, for many use- 
ful suggestions which improved the paper substantially. We thank Mike Vogeley, Ann 
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and by the Smithsonian Institution. R.M. acknowledges NASA grant NGT-50819. 
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TABLE 1 



PREVIOUS MEASUREMENTS OF (712 





oaiiipie 


AT , 


(7i2l^KmS j 


Geller & Peebles 1973 


RCl 


226 


206 


Peebles 1976 


RCl 


422 


200" 


Peebles 1979 


KOS^ 


166 


500 


Davis & Peebles 1983 


cfAr 


2397 


340±40 


Roan ci ril 1 OlU'X 

Dean ex at. iyoo 


A A R Q<i 
AArtD 


ozu 


ZOUit OU 


Hale-Sutton et al. 1989 


SAAO^ 


264 


600 ± 140 


Mo et al. 1993 


Several 


~ 700 -2400 


~ 200 -1000 


Fisher et al. 1994 


IRAS 1.2 Jy 


5313 


317 ±45 


Zurek et al. 1994 


CfAl 


2397 


~ 300-600 


'^"^Excludes the Virgo 
'^^^Kirshner, Oemler & 


cluster. Including Virj 
Schechter 1978 


?o, (712 ~ 600 km 





(^)Huchra et al. 1983 



Peterson et al. 1983 
(^)Metcalfe et al. 1989 
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TABLE 2 



SURVEY SAMPLES 



Sample 




il (ster) 


Volume (10^ h-^Mpc^) 


CfA2 North 


6480 


1.23 


7.1 


CfA2 South 


2741 


0.77 


4.4 


SSRS2 


3591 


1.12 


6.5 


CfA2+SSRS2 


12812 


3.12 


18.0 
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TABLE 3 

CORRELATION FUNCTION ^(r) 

Sample tq (/j""*" Mpc) 7 x^/d.o.f. 



CfA2 North 
No Clusters 



5.83±0.18 
4.77±0.17 



-1.80±0.03 
-1.80±0.04 



11.5/18 
16.8/18 



CfA2 South {b > W) 
No Clusters 



4.75±0.20 
3.74±0.18 



-1.99±0.05 
-2.11±0.08 



18.9/18 
12.0/18 



SSRS2 
No Clusters 



5.08±0.23 
4.77±0.21 



-1.85±0.06 
-1.80±0.05 



15.1/18 
12.2/18 



CfA2+SSRS2 
No Clusters 



5.97±0.15 
4.95±0.13 



-1.81±0.02 
-1.73±0.03 



19.8/18 
26.4/18 
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TABLE 4 



PAIRWISE DISPERSION FOR CfA2 NORTH 



Clommpnts 


r-p 




Y^/dof 







Y^/dof 







{h-^ Mpc) 


(km s"-*-) 






(km s"-*-) 








0.15 


190±27 


9.7/13 


0.72 


578±77 


7.8/13 


0.86 




0.3 


259±27 


5.2/13 


0.97 


592±88 


4.6/13 


0.98 


No Infall 


0.6 


305±37 


8.5/13 


0.81 


770±116 


4.9/13 


0.98 




1.2 


274±26 


10.1/13 


0.69 


697±92 


5.0/13 


0.98 




2.4 


226±29 


7.3/13 


0.88 


510±85 


6.6/13 


0.92 




4.8 


79±37 


7.4/13 


0.88 


149±89 


7.5/13 


0.88 




0.15 


199±27 


10.2/13 


0.68 


509±75 


7.8/13 


0.86 




0.3 


270±26 


5.3/13 


0.97 


640±80 


4.6/13 


0.98 


Infall 


0.6 


329±35 


8.8/13 


0.79 


825±124 


4.8/13 


0.98 




1.2 


317±26 


11.2/13 


0.60 


840±103 


4.6/13 


0.98 




2.4 


307±26 


8.4/13 


0.82 


751±93 


6.9/13 


0.91 




4.8 


238±19 


5.5/13 


0.96 


588±61 


5.2/13 


0.97 




0.15 


182±27 


2.2/13 


1.00 


384±82 


2.4/13 


1.00 




0.3 


128±22 


20.0/13 


0.09 


222±55 


22.3/13 


0.05 


No Infall 


0.6 


157±26 


4.6/13 


0.98 


333±74 


6.1/13 


0.94 


No Clusters 


1.2 


128±21 


6.6/13 


0.92 


276±57 


7.3/13 


0.89 




2.4 


109±27 


4.4/13 


0.99 


219±68 


4.7/13 


0.98 






7-1-1 ^fi 




u .yy 






n 1 9 




0.15 


189±27 


2.3/13 


1.00 


407±104 


2.4/13 


1.00 




0.3 


143±20 


16.3/13 


0.23 


258±48 


18.8/13 


0.13 


Infall 


0.6 


180±24 


4.1/13 


0.99 


392±70 


5.6/13 


0.96 


No Clusters 


1.2 


170±18 


6.6/13 


0.92 


388±56 


7.3/13 


0.89 




2.4 


191±19 


3.0/13 


1.00 


435±61 


3.5/13 


1.00 




4.8 


197±20 


8.9/13 


0.78 


452±61 


12.6/13 


0.48 



(") p(y«2) from Equation 18, assuming isotropic exponential P{vi). 
'^^^ Isotropic exponential Pivl^) (Equation 15). 
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TABLE 5 



PAIRWISE DISPERSION FOR CfA2 SOUTH 



Clommpnts 


r-p 

{h-^ Mpc) 


(km s"-*-) 


Y^/dof 





(km s"-*-) 


Y^/dof 







0.15 


166±26 


6.0/13 


0.95 


400±78 


4.9/13 


0.98 


No Infall 


0.3 


198±35 


5.0/13 


0.98 


457±119 


3.7/13 


0.99 




0.6 


109±15 


16.3/13 


0.23 


244±45 


15.0/13 


0.31 




1.2 


147±23 


8.9/13 


0.78 


348±69 


5.0/13 


0.97 




2.4 


109±29 


5.6/13 


0.96 


243±76 


4.3/13 


0.99 




4.8 


63±85 


5.6/13 


0.96 


115±197 


5.7/13 


0.96 




0.15 


179±28 


6.3/13 


0.93 


366±80 


4.7/13 


0.98 




0.3 


207±34 


5.0/13 


0.97 


461±104 


3.6/13 


1.00 


Infall 


0.6 


131±14 


16.9/13 


0.21 


295±45 


14.4/13 


0.34 




1.2 


197±23 


10.9/13 


0.62 


463±72 


4.6/13 


0.98 




2.4 


199±24 


8.9/13 


0.78 


483±80 


3.7/13 


0.99 




4.8 


198±35 


3.1/13 


1.00 


475±111 


4.8/13 


0.98 




0.15 


130±22 


0.4/8 


1.00 


254±65 


1.0/8 


1.00 




0.3 


164±39 


4.1/8 


0.85 


353±111 


3.4/8 


0.91 


No Infall 


0.6 


92±19 


1.9/8 


0.98 


185±53 


2.1/8 


0.98 


No Clusters 


1.2 


78±24 


0.4/8 


1.00 


162±61 


0.4/8 


1.00 




2.4 


15±210 


2.6/8 


0.96 


27±472 


2.7/8 


0.95 








m 


/ 






m 








0.15 


135±21 


1.1/8 


0.99 


264±52 


1.0/8 


1.00 




0.3 


170±40 


4.1/8 


0.85 


160±118 


11.6/8 


0.17 


Infall 


0.6 


110±17 


2.0/8 


0.98 


225±48 


2.0/8 


0.98 


No Clusters 


1.2 


118±18 


1.0/8 


1.00 


253±54 


0.4/8 


1.00 




2.4 


108±22 


4.7/8 


0.78 


256±70 


3.6/8 


0.89 




4.8 


150±41 


1.8/8 


0.99 


330±113 


2.1/8 


0.98 



(") p(y«2) from Equation 18, assuming isotropic exponential P{v\). 
'^^^ Isotropic exponential P{vyi) (Equation 15). 
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TABLE 6 



PAIRWISE DISPERSION FOR SSRS2 



Clommpnts 


r-p 




Y^/dof 







Y^/dof 







{h-^ Mpc) 


(km s"-*-) 






(km s"-*-) 








0.15 


150±32 


3.4/13 


1.00 


329±70 


3.3/13 


1.00 


No Infall 


0.3 


105±30 


10.2/13 


0.68 


199±86 


11.3/13 


0.59 




0.6 


131±23 


10.2/13 


0.68 


287±66 


5.1/13 


0.97 




1.2 


109±28 


4.7/13 


0.98 


210±73 


5.8/13 


0.95 




2.4 


48±45 


3.8/13 


0.99 


95±98 


3.8/13 


0.99 




4.8 


11±151 


13.0/13 


0.45 


9±484 


20.9/13 


0.08 




0.15 


160±29 


3.6/13 


0.99 


334±95 


3.3/13 


1.00 




0.3 


115±26 


8.2/13 


0.83 


210±69 


9.2/13 


0.76 


Infall 


0.6 


155±21 


5.4/13 


0.96 


345±65 


5.2/13 


0.97 




1.2 


152±22 


3.8/13 


0.99 


324±65 


5.5/13 


0.96 




2.4 


156±21 


7.8/13 


0.86 


367±62 


5.2/13 


0.97 




4.8 


189±28 


9.0/13 


0.78 


453±80 


7.1/13 


0.89 




0.15 


120±32 


3.1/8 


0.93 


227±72 


3.7/8 


0.88 




0.3 


95±25 


3.9/8 


0.86 


179±67 


4.6/8 


0.80 


No Infall 


0.6 


109±23 


3.8/8 


0.87 


226±65 


4.5/8 


0.81 


No Clusters 


1.2 


82±33 


2.2/8 


0.97 


161±65 


2.3/8 


0.97 




2.4 


43±213 


3.4/8 


0.91 


85±59 


2.3/8 


0.97 
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o.Of o 


n S7 

U .O I 




^A) j O 


n sn 
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0.15 


131±27 


2.8/8 


0.94 


241±64 


3.4/8 


0.91 




0.3 


109±22 


3.1/8 


0.93 


213±61 


3.7/8 


0.88 


Infall 


0.6 


131±20 


3.9/8 


0.87 


283±63 


4.5/8 


0.81 


No Clusters 


1.2 


133±23 


1.8/8 


0.99 


275±64 


1.9/8 


0.98 




2.4 


142±20 


5.9/8 


0.66 


317±59 


4.4/8 


0.82 




4.8 


126±23 


5.0/8 


0.75 


277±64 


4.9/8 


0.77 



(") p(y«2) from Equation 18, assuming isotropic exponential P{v\). 
'^^^ Isotropic exponential P{vyi) (Equation 15). 
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TABLE 7 



PAIRWISE DISPERSION FOR CfA2 + SSRS2 



L^l i 11 i iCli u o 


r-p 


at 


/Hof 


o 


^2 


/Hof 


o 




{h-^ Mpc) 


(km s"-*-) 






(km s"-*-) 








0.15 


233±22 


4.9/13 


0.98 


545±65 


1.4/13 


1.00 


No Infall 


0.3 


224±21 


9.2/13 


0.76 


540±57 


5.6/13 


0.96 




0.6 


174±16 


28.5/13 


0.01 


534±68 


16.1/13 


0.24 




1.2 


209±18 


13.8/13 


0.39 


515±55 


5.7/13 


0.96 




2.4 


172±20 


14.0/13 


0.38 


398±54 


8.5/13 


0.81 




4.8 


9±231 


14.0/13 


0.38 


111±93 


4.7/13 


0.98 




0.15 


239±19 


5.1/13 


0.97 


564±65 


1.5/13 


1.00 




0.3 


238±21 


9.8/13 


0.71 


579±71 


5.6/13 


0.96 


Infall 


0.6 


200±16 


32.7/13 


0.01 


599±64 


15.6/13 


0.27 




1.2 


260±17 


17.0/13 


0.20 


647±57 


5.0/13 


0.97 




2.4 


274±17 


20.2/13 


0.09 


654±56 


6.4/13 


0.93 




4.8 


248±16 


7.8/13 


0.86 


601±51 


4.6/13 


0.98 




0.15 


168±26 


3.0/13 


1.00 


345±77 


3.3/13 


1.00 




0.3 


139±19 


16.4/13 


0.23 


269±48 


19.0/13 


0.12 


No Infall 


0.6 


128±17 


8.2/13 


0.83 


270±48 


9.0/13 


0.78 


No Clusters 


1.2 


109±14 


9.1/13 


0.77 


230±39 


9.7/13 


0.72 




2.4 


75±21 


5.0/13 


0.98 


156±50 


4.9/13 


0.98 




A S 


Q-i-/ini 


o.u/ io 


U.yo 


1 Q_l_ 

ioitozo 


DZ.Z/ io 


u.uu 




0.15 


176±23 


3.2/13 


0.99 


380±62 


0.1/13 


1.00 




0.3 


158±17 


13.2/13 


0.43 


304±43 


16.3/13 


0.23 


Infall 


0.6 


156±15 


8.1/13 


0.83 


337±45 


8.2/13 


0.83 


No Clusters 


1.2 


153±12 


8.5/13 


0.81 


351±38 


8.1/13 


0.83 




2.4 


181±13 


7.8/13 


0.86 


418±39 


3.4/13 


1.00 




4.8 


190±19 


4.4/13 


0.99 


444±56 


6.5/13 


0.93 



(") p(y«2) from Equation 18, assuming isotropic exponential P{vi). 
'^^^ Isotropic exponential Pivl^) (Equation 15). 
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TABLE 8 



BEST ESTIMATES OF (712 



Sample cri2(< 0.8) [kms ^] Formal Error [kms ^] External Error [kms ^] 

CfA2 North 647 52 

CfA2 South 367 38 

SSRS2 272 42 

CfA2+SSRS2 540 36 180 



38 




Figure 1: Pairwise velocity distribution functions from Equations 15 and 18. 
Curves are labeled by the equation number. Equation 15 describes the exponen- 
tial pairwise distribution. Equation 18 describes the pairwise velocity distribu- 
tion generated by an isotropic exponential distribution of single-galaxy velocities. 



39 



a) Cooled Veloolty Field 
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c) Isotropic Exponential P(vi) 
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Figure 2: C(rp,7r) for the N-body simulation. Solid contour indicates ^ = 1. 
Countour intervals are = 0.1 for ^ < 1 and Alog^ = 0.1 for ^ > 1. Dotted 
contours indicate ^ < 0. (a) particle positions and velocities taken directly from 
the cooled simulation (see text), (b) pure Hubble flow - peculiar velocities set 
to zero, (c) particle velocities replaced with random deviates drawn from the 
isotropic exponential distribution in Equation 16, (d) best-fltting model from 
Equation 8. 
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Figure 3: The projected CF, w{rp), determined from numerical integration of 
^(rpjir) out to TTmax- The solid line is the input model: tq = 5.8 /i""*" Mpc and 
7 = —1.8. C(rp,7r) comes from Equation 8. Squares are for iTmax — QOh~^ Mpc, 
filled circles are for iTmax — 30 /i""*" Mpc. For the data, we integrate to iTmax — 
30 /i-^ Mpc. 
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Figure 4: Projected correlation function w{rp) for the N-body simulation. The 
solid line is the power-law fit. 
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Figure 5: Two-point correlation function from pair counts in the N-body sim- 
ulation. Solid line is the inversion of the power-law fit to w{rp). 
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Figure 6: Distribution of pairwise velocities from Monte-Carlo simulations. Par- 
ticle positions are from the N-body simulation. N-body velocities are replaced 
with velocities drawn randomly from the distribution of single-galaxy velocities 
in Equation 15. Magnitudes of the velocities are exponentially distributed with 
(Ji = 250kms~^; directions are random. The solid line is the distribution of 
pairwise velocities predicted by Equation 18. The dashed line is the exponential 
distribution of 1-d pairwise velocities with the same second moment, (J12 = 2(Ji. 
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Figure 7: Fits to ^{rp, it) derived from N-body simulation with velocities drawn 
from an isotropic exponential distribution (see figure 6). Numbers in the upper 
right-hand corner indicate the range of Vp in /i~^Mpc. Note that the vertical 
scale is dilferent for each range of Vp . 
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Figure 8: Solid lines are the first two central moments of the distribution P{vi2) 
for the N-body simulation. Open squares represent the fitted (J12 ignoring infall. 
The dotted line is the measured (J12 using the DP83 infall model (Equation 19) 
with F — 1.5, which is shown here as the dashed line. 
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Figure 9: ^{rp, it) for the four samples in Table 2 (MV weighting). The contour 
levels are described in the caption for Figure 2. 
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Figure 10: C(rp,7r) for the four samples in Table 2 (uniform weighting). The 
contour levels are described in the caption for Figure 2. 
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Figure 11: w{rp) for the four samples in Table 2. The solid line is the power-law 
fit given in Table 3. 
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Figure 12: Fits to ^{ir) at different Vp for CfA2N. Solid line is the pairwise 
exponential; dashed line represents the single-galaxy exponential. Numbers in 
the upper right-hand corner indicate the range of Vp in Mpc. These fits do 
not include a correction for infall. 
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Figure 13: Fits to ^{ir) at different Vp for CfA2S. Solid line is the pairwise 
exponential; dashed line represents the single-galaxy exponential. Numbers in 
the upper right-hand corner indicate the range of Vp in Mpc. These fits do 
not include a correction for infall. 
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Figure 14: Fits to ^{ir) at different Vp for SSRS2. Solid line is the pairwise 
exponential; dashed line represents the single-galaxy exponential. Numbers in 
the upper right-hand corner indicate the range of Vp in Mpc. These fits do 
not include a correction for infall. 
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Figure 15: Fits to ^{ir) at different Vp for CfA2+SSRS2. Solid line is the pairwise 
exponential; dashed line represents the single-galaxy exponential. Numbers in 
the upper right-hand corner indicate the range of Vp in Mpc. These fits do 
not include a correction for infall. 
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Figure 16: Scaling of (J12 with Vp. The samples are CfA2N (open squares), 
CfA2S (filled squares), SSRS2 (open circles) and CfA2+SSRS2 (filled circles). 
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Figure 17: Fits of Equations 15 and 18 to C(rp,7r) for Vp < l/i~^Mpc in 
CfA2+SSRS2. Solid line is the pairwise exponential; dashed line represents the 
single-galaxy exponential. 
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Figure 18: ^{rp, it) for the four samples in Table 2 after removing Abell clusters 
with R > 1. The contour levels are described in the caption for Figure 2. 
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Figure 19: Comparison with (J12 derived from N-body simulations. The filled 
squares are our results for CfA2+SSRS2 with no correction for infall. Open 
squares represent CfA2+SSRS2 with Abell clusters removed. Filled circles are 
from Davis et al (1985), where we have taken the mean of the two simula- 
tions from their figure 15. CC comes from Couchman & Carlberg (1992). CO 
represents Cen & Ostriker (1993); CO represents their "bright" sample, which 
includes only galaxies with mass greater than 1.1 x 10^° Mq. 




Figure 20: Cumulative distribution of one-dimensional velocity dispersions 
n{> a) for systems of galaxies. The curves labeled by erg are the predictions of the 
Press-Schechter prescription; erg is the normalization of a flat CDM power spec- 
trum. Curves labeled by a represent the range of distributions 
n{> a) oc lO"'^ derived from optical surveys. X-ray observations correspond 
best with erg = 0.6 for flat CDM. All distributions are normalized to match the 
observed number density of Abell clusters with richness class R>1. 
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Figure 21: Convergence of (J12 averaged over virialized systems. Umax is the 
upper limit of integration in Equation 22. Curves correspond to the distributions 
n{> a) in figure 20 and are labeled accordingly. 
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Figure 22: Monte-Carlo analysis of the error in (J12 as a function of survey vol- 
ume. The pairwise velocity dispersion is biased low when the volume does not 
include a significant number of clusters with a > Uconv, where Uconv is the ve- 
locity dispersion a at which (J12 converges (see figure 21). The dispersion Umax 
is the upper limit of integration in Equation 22. The slopes a of the exponential 
n{a) correspond to the observed distributions in figure 20. Arrows indicate (from 
left to right) the volumes of CfAl, CfA2N and CfA2+SSRS2. 
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